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OBJECT 


The  object  of  the  effort  covered  in  this  report  is  to  analyze 
the  stresses  induced  in  the  plastic  rotating  band  during  the  time  in 
which  the  projectile  travels  inside  the  gun  after  it  has  been  en- 
graved. In  this  period  the  acceleration  of  the  projectile  causes 
radial  and  circumferential,  as  well  as  axial  normal  stresses.  Shear- 
ing stresses  are  also  generated  in  two  directions.  The  assessment  of 
the  levels  of  stresses  in  the  rotating  band  not  only  will  furnish 
necessary  information  in  the  selection  of  the  plastic  materials,  it 
also  will  shed  light  on  the  gun  wear  problem  because  the  contact 
pressure  between  the  surfaces  of  the  gun  and  the  rotating  band  is 
also  determined  by  the  analyses. 

An  axisymmetric  model  of  the  rotating  band  in  flight  inside  the 
gun  after  it  has  passed  the  forcing  cone  has  been  formulated  and 
analyzed  mathematically.  The  analytical  solutions  of  stress  compo- 
nents were  then  used  in  a sample  calculation  based  on  a sample 
material  with  E = 6895  MPa  (1,000,000  psi)  and  a sample  geometry  of 
a 105  mm  gun.  Results  show  that  the  maximum  compressive  stress  at 
the  interface  of  the  gun  and  the  rotating  band  is  about  110  MPa 
(16,000  psi)  whereas  at  the  band  seat  a tensile  stress  of  about  the 
same  magnitude  is  induced.  Shearing  stresses  are  only  of  the  order 
of  0.7  MPa  (several  hundred  psi).  It  is  observed  that  due  to  the 
assumption  of  the  rigidity  of  the  gun  and  the  shell  body  the  normal 
stresses  are  higher  than  the  actual  stresses.  The  stress-free  con- 
ditions assumed  for  the  circumferential  and  axial  displacements 
cause  the  shearing  stress  to  be  too  low.  Refinements  of  the  mathe- 
matical model  currently  in  use  will  not  be  a major  undertaking  and 
are  recommended. 


INTRODUCTION 


In  the  design  of  a plastic  rotating  band,  the  design  analysis 
can  be  arbitrarily  divided  into  three  separate  stages.  The  first 
stage  is  to  deal  with  the  engraving  process.  As  the  rotating  band 
enters  the  forcing  cone  the  outside  diameter  of  the  band  is  com- 
pressed at  first  and  then  engraved  by  the  lands  on  the  gun  tube. 

The  objective  of  the  analysis  of  this  stage  of  action  is  to  deter- 
mine the  pressure  exerted  on  the  gun  tube  by  the  rotating  band  as  it 
is  being  engraved.  The  second  stage  of  analysis  deals  with  the  dyna- 
mics and  the  stress  problems  occurring  after  the  projectile  has 
emerged  from  the  forcing  cone  and  during  its  travel  through  the  gun 
tube.  It  ends  at  the  instant  the  projectile  leaves  the  muzzle  of 
the  gun.  The  third  stage  refers  to  the  problems  relating  to  the  band 
during  the  free  flight  of  the  projectile  (Fig  1). 

After  the  rotating  band  has  been  engraved  by  the  forcing  cone 
the  projectile  is  made  to  accelerate  through  the  gun  both  axially 
and  circumferentially,  acquiring  the  desired  muzzle  velocity  and 
spin  at  its  completion  of  the  travel  inside  the  gun.  During  this 
period,  pressure  is  generated  between  the  mating  surfaces  of  the 
band  and  the  barrel  due  to  the  inertial  forces  acting  on  the  band. 

This  pressure  has  direct  influence  on  the  frictional  force  which 
exists  at  the  interface  of  the  radial  contact.  It  therefore  deter- 
mines the  wear  of  the  mating  surfaces.  Besides  the  radial  contacts 
there  are  the  lateral  contacts  between  the  band  and  the  gun  barrel. 
Very  high  contact  pressure  exists  there  due  to  the  angular  accelera- 
tion of  the  projectile.  These  pressures  and  the  stresses  induced  in 
the  band  during  the  post  engraving  period  are  the  subject  of  analysis 
of  this  paper. 

The  actual  problem  of  stress  analysis  of  the  rotating  band  dur- 
ing its  travel  in  the  gun  is  a complicated  one.  Foremost  among  the 
complexities  involved  is  the  geometry  of  the  mating  surfaces  of  the 
band  in  contact  with  the  gun.  These  surfaces  are  in  the  form  of 
lands  and  grooves  after  the  band  has  been  engraved.  Since  the  gun 
barrel  and  the  shell  body  both  deform,  though  to  a much  lesser  ex- 
tent compared  to  the  deformation  of  the  band,  to  analyze  the  stresses 
in  the  band  one  needs  also  to  analyze  the  deformations  of  the  gun  and 
the  shell.  The  complete  formulation  thus  involves  three  problems 
(the  gun,  the  band,  and  the  shell)  coupled  together.  Simplifications 
must  be  made  before  analysis  can  proceed. 

In  the  following  analysis  the  rotating  band  is  considered  to  be 
a smooth  cylindrical  shell  instead  of  the  serrated  configuration. 

This  assumption  ignores  the  lateral  surfaces  of  the  land  and  thus 
the  stress  variations  due  to  the  irregularities  in  the  circumferential 
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direction.  Such  calculation  should  give  a meaningful  estimate  of  the 
order  of  magnitude  of  the  stress  levels  at  the  contact.  This  estima- 
tion of  stress  level  would  be  useful  in  the  selection  of  materials 
suitable  for  the  application  with  regard  to  their  strengths  and  wear 
resistances.  A second  assumption  is  to  ignore  the  deformations  of 
the  gun  and  the  shell  so  as  to  decouple  the  problem  of  the  band  de- 
formation. This  means  that  in  the  analysis  we  either  take  the  de- 
formations of  the  gun  and  the  shell  to  the  zero  or  use  some  predeter- 
mined values  for  these  deformations,  say  from  some  experiments. 

Besides  the  assumptions  made  above,  it  will  be  assumed  that  the 
rotating  band  under  the  given  load  is  stressed  within  the  elastic 
limit.  Therefore  the  subsequent  analysis  will  be  based  on  the  theory 
of  elasticity. 


MATHEMATICAL  ANALYSIS 


According  to  the  theory  of  elasticity  three  equations  of  motion 
in  the  cylindrical  coordinates  r,  0,  and  z are  as  follows  (Ref  1): 
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and  u,  v,  w are  the  displacements  in  the  radial  (r) , circumferential 
(0) , and  axial  (z)  directions  of  a point  on  the  band,  referring  to  a 
frame  of  reference  rotating  with  the  band,  and  ar,  aQ,  az  are  the 
acceleration  components. 
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Let  the  angular  velocity  and  the  angular  acceleration  be  to  and  <o, 
the  rifling  angle  be  6 and  the  outer  radius  of  the  band  be  b.  The 
components  of  the  acceleration  of  a point  on  the  band  as  observed 
from  a laboratory  frame  of  reference  are  calculated  as: 


ar  - -r„2-2„  £ ♦ 2-5 


n2 
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To  further  facilitate  the  analysis,  the  band  is  assumed  to  be  in 
plane  strain;  thus  the  z coordinate  can  be  ignored.  It  is  further 
observed  that  the  displacements  are  independent  of  their  angular 
positions.  This  means  that  the  0 coordinate  can  also  be  ignored. 

The  remaining  independent  variables  are  the  radial  distance  r and 
the  time  t. 

Using  all  the  above  assumptions  and  substituting  Equation  (3) 
into  Equation  (1)  we  obtain  the  one-dimensional  time-dependent  model 
of  the  rotating  band  as  described  by  the  following  three  partial 
differential  equations. 
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It  can  be  seen  that  the  first  and  the  second  equations  of 
Equation  (4)  are  coupled  through  the  Coriolis  terms. 


Corresponding  to  Equation  (4)  a system  of  boundary  conditions 
must  be  prescribed  to  depict  the  physical  conditions  existing  at  the 
inner  and  the  outer  contact  surfaces  of  the  band.  For  the  radial 
displacement  the  boundary  conditions  state  that  its  values  must 
vanish  or  be  equal  to  some  predetermined  constants.  These  constants 
can  be  selected  from  some  known  data  pertaining  to  the  particular 
ordnance  system  in  consideration,  i.e.,  u(a)  = u(b)  = 0,  or  u(a)  = A, 
u(b)  = B,  where  a and  b are  the  inner  and  the  outer  radii  of  the  band, 
and  A and  B are  the  predetermined  values  of  deformation  of  the  inner 
and  the  outer  surfaces. 

For  the  circumferential  displacement  v the  boundary  conditions 
. are  v(a)  = 0,  which  corresponds  to  the  condition  of  no  relative 

motion  at  the  interface  between  the  band  and  the  shell,  and  v(b) 
being  unspecifiable.  At  the  contact  between  the  band  and  the  gun 
either  the  interface  can  be  assumed  to  be  stress  free  or  a predeter- 
mined frictional  stress  can  be  assumed.  The  boundary  conditions  for 
the  axial  components  w are  similar  to  those  of  v. 

t 

In  Equation  (4)  the  angular  velocity  and  the  angular  accelera- 
tion are  given  functions  of  time.  For  each  ammunition  system,  de- 
pending on  the  condition  of  the  gun  and  the  zone  of  firing,  the 
interior  ballistic  performances  vary.  Some  typical  data  will  be  used 
, in  the  analysis  to  get  representative  answers  which  will  disclose  the 

stress  levels  induced  in  the  particular  system. 

The  individual  equations  in  Equation  (4)  are  similar  in  their 
structure.  To  determine  the  displacements  due  to  the  imposed  in- 
ertial forces,  one  must  solve  for  the  particular  solutions  due  to  the 
inertial  forces.  To  accomplish  this  the  homogeneous  solutions  must 
be  obtained  first.  This  can  be  accomplished  by  the  separation  of 
variables  technique. 
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where 


e2  - 


A + 2y 


p a2 


The  homogeneous  part  of  Equation  (6)  is  the  Bessel's  equation 
of  the  first  order,  whose  solution  can  be  written  in  the  form 


U(r)  = Aj  Jj(Br)  + ^ Yj(Br)  (7) 

where  Jj  and  Yj  are  the  Bessel's  functions  of  the  first  and  the 
second  kind.  The  subscripts  denote  the  order  of  the  Bessel's  func- 
tion. 


Upon  imposing  the  boundary  conditions 


U(a)  - U(b)  = 0 


(8) 


the  eigenvalue  equation  is  obtained. 


Jj(Ba)  Yj(Bb)  - JjCBb)  Yj(Ba)  = 0 (9) 


If  a new  variable  B = Ba  and  a parameter  k = b/a  are  introduced, 
then  Equation  (9)  can  be  cast  in  the  form  of 

J (B)  Y (kB)  - J (k?)  Y (B)  = 0 (10) 


Since  there  is  no  available  information  regarding  the  eigen- 
values of  this  equation,  numerical  solutions  will  be  performed  to 
obtain  the  eigenvalues.  These  values  are  needed  for  subsequent 
calculations . 


Assuming  the  eigenvalue^  are  Bn,  the  set  of  orthogonal  functions 
corresponding  to  the  set  of  Bn  are  given  by 


= Jl<*n  £>  - K V6n 


where 


K 


a relation  derivable  from  Equation  (10). 
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To  solve  the  inhomogeneous  problem.  Green's  function  method  will 
be  employed.  Let  <j>(r,t)  be  a general  representation  of  the  inhomo- 
geneous terms  in  Equation  (4).  The  particular  solution  in  terms  of 
4>  can  be  represented  by  the  following  integral. 


u(r,t) 


t rb 


G(r,C,t-T)<J>(£,T)d£dT 


0 a 


where  G(r,£,t-T)  is  Green's  function  of  the  partial  differential 
equation. 

Substituting  Equation  (12)  in  the  first  equation  of  Equation  (4) 
yields 


(A  + 2p)  | | 


<KS,T)d£dT  = <Hr,t) 


from  which  it  can  be  seen  that 

(X+2p)  -p^f  = «(r-C)«(t-T)  (13) 

3 1 

in  which  the  right  hand  side  is  the  product  of  two  Dirac-delta 
functions. 

To  determine  Green's  function  G(r.C.t-t)  it  is  necessary  to  ex- 
pand the  delta  function  in  the  spatial  variable  into  an  infinite 
series  in  the  orthogonal  function  ifi(8  *1) . 


Letting 


(r-C)  = l A (OSlPCB 
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and  due  to  the  orthogonality  relation  existing  between  i/j- functions 


corresponding  to  different  Bn>  i.e.. 
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we  obtain 
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To  separate  variables  in  Equation  (13)  we  assume  that 

G (r , C > t-  t)  = l ^Bn(t-x)5*(Bn  f)iKBn  f)  (16) 
n=l 


Substituting  Equation  (16)  into  Equation  (13)  and  using  the  fact 
that  (tKJ  — ) satisfies  the  homogeneous  part  of  Equation  (6)  yields  an 
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Equation  (17)  has  the  following  solution 
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From  some  data  available  for  the  105  mm  system  a typical  angular 
velocity  w vs  t relationship  can  be  represented  by  a straight  line, 
namely. 


uj(t) 


(21) 


where  ioQ  is  the  angular  velocity  of  the  shell  at  the  exit  time  tQ. 

Using  the  expression  of  gj(t)  in  Equation  (21)  and  performing 
the  integration  with  respect  to  the  t-variable,  the  time  part  of  the 
integral  in  Equation  (20)  is  obtained. 
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The  integration  with  respect  to  the  variable  % is  less  straight- 
forward and  is  therefore  worked  out  and  presented  with  a few  details 
in  the  Appendix.  Substituting  Equations  (22),  (A4) , and  (A6)  into 
Euqation  (20)  yields 


u(r , t)  = l Cn^(6nr)T  (aRt)  (23) 

n=  1 


where 


2 oo. 
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with 
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In  Equation  (25)  K = b/a  and  T*  is  the  function  inside  the  bracket  of 
Equation  (22) , 
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The  Circumferential  Displacement,  v(r,t) 

Let  v(r,t)  = V(r)sinyt  (26) 


Substituting  Equation  (26)  into  the  second  equation  of  (4),  with 
the  Coriolis  term  dropped,  yields 

^7  * 7 HK  * ^2'7)V  - 0 (27) 

2 2 

where  k = — — . 
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The  boundary  condition  at  the  interface  between  the  rotating 
band  and  the  shell  body  requires  that  the  circumferential  displace- 
ment v there  vanishes  to  correspond  to  the  condition  of  positive 
anchorage  of  the  band  to  the  shell.  At  the  outer  radius  r = b it  is 
assumed  that  the  shearing  stress  is  zero  to  correspond  to  a smooth 
contact  without  friction.  However,  since  the  frictional  stress  is 
dependent  on  the  radial  pressure  at  the  interface,  to  introduce  fric- 
tion means  to  introduce  a coupling  between  the  radial  and  circum- 
ferential displacements  through  the  boundary  condition.  This  would 
make  the  mathematical  problem  unnecessarily  complicated  for  the  pur- 
pose of  this  analysis.  A compromise  approach  would  be  to  use  a pre- 
determined value  of  friction  and  thus  introduce  an  inhomogeneous 
boundary  condition.  In  this  paper  the  zero  stress  condition  will  be 
used  for  mathematical  expediency  with  the  anticipation  that  the  re- 
sulting shear  stress  may  thus  be  lower  than  the  actual  value. 


The  solution  of  Equation  (27)  is 


V (r)  = Bj  J1(<r)+B2 


(29) 


From  the  first  boundary  condition  in 


Equation  (28)  it  follows  that 


B1  J1C<a)+B2  Y1(<a)  = 0 


or 


B-  " ' Y‘{1eaT  B1  * 


(30) 


The  second  boundary  condition  leads  to  the  following 

Bx[bK  J1’(<b)-J1(<b)]+B2(bic  Y0'(Kb)-Y0(Kb)]  = 0 (31) 

Lgenvalue  equation,  after  simplification,  becomes 


Thus,  the  eij 


J ' (<b)Y1(Kb)-Y1' (<b)J1(<b)  = 0 . 
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Equation  (32)  can 
the  derivatives  of  the 


be  further  reduced  by  using 
Bessel's  functions. 


the  identities  of 
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then  Equation  (32)  becomes 
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This  equation  -ill  be  solved  numerically  to  determine  the  values  of 
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The  eigenfunctions  are  given  by 
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Using  <f>(r,t)  = prw  in  evaluating  the  particular  integral  v,  we  have 
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The  integration  of  the  time  part  gives 


T (Y-t)  = (1-cosy  t) 

n YnZ0  n 


The  integration  with  respect  to  £ yields,  after  some  elementary  com- 
putation, the  following  expression  for  v. 


v(r>  t)  - l D <H«.  r)T*(«  t) 
n=l  n n n 
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In  Equation  (38)  k = b/a  and  T*  is  the  function  inside  the  paren- 
theses of  Equation  (35), 

*0<Kna)  ” J0(Kna)‘MnYOf,:na)- 

" Jl<’tnb>-M„Yl''<nb>- 
^2(Kna)  = J2(Kna)-MnY2(<na) , etc. 
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The  Axial  Displacement,  w(r,t) 

Let  w(r,t)  = W(r)  sin  pt  (3?T 

Substituting  Equation  (39)  into  the  third  equation  of  Equation  (4) 
yields 

~7  + F 37  + (q2)W  = 0 (40) 


where 


.2  _ pp 


The  considerations  for  the  boundary  conditions  are  identical 
with  those  for  the  circumferential  displacement  v.  Therefore, 

W(a)  = 0 (4] 


37(b)  = 0. 

The  solution  of  Equation  (40)  is 


W(r)  = E1J0(qr)+E2Y0(qr) 


From  the  first  boundary  condition  in  Equation  (41)  it  follows  that 


ElJ0(cla)+E2Yo(qa^  = 0 


E - J°(qa)  p 

2 ^(qa)  Ei  • 


The  second  boundary  condition  leads  to  the  following. 


Vo  (qb)+E2Y0,(qb)  = 0 


13 


The  eigenvalue  equation  in  this  case  is 


J0  ^b)Y0^a)-Jo(qa)Y0,(qb)  = 0 


and  the  eigenfunctions  are 


!>  ' V%  ?'-Ln  V«n 


where 


L . W> 

~ "t 7~ T^T 


n J0Cqna)  ' 


With  0 = P ( 0) , we  have 


f*  rb  » ?<Kq  |)<Hq  §)  b . 

- J Jx  p-fTp^ (ptt^)sin  Pn(t-T)d5dx 


The  integration  of  the  time  part  gives 


T<V>  ■ Ttpyt^  C1-«SV)  (4' 

The  integration  with  respect  _to  £ yields,  after  some  calculation, 
the  following  expression  for  w. 


, b)-^  (q  a) 

nk,q_)  = ~ vV  ^ (51 

k^0Z(qnb)-^12(qna) 

where  ^ (^b)  = JQ  (qnb)  - (q^)  , 

^l(qna)  = J1Cqna)-LnYo(qna).  etc. 

Stress  Components 

The  rotating  band  is  assumed  to  undergo  elastic  deformation; 
therefore,  the  stress-strain  relation  is  given  by 


Tkl  ~ Ae6k£+2yek£  (51) 

According  to  the  pattern  of  deformation  assumed  in  this  analysis  the 
components  of  strain  are 


err 

9u 

9r  » 

er0 

1 r9v 
2l9r 

1 

eee 

u 

r ’ 

erz  = 

1 9w 

2 9r 

e = 
zz 

0 , 

e0  z = 

0 

1 

The  corresponding 

stress  components 

are: 

• 

Trr 

(A+2y)|H  + 

4 

* Tr0 

Tee 

(X+2p)^  + 

. 9u 
9r 

’ Trz 

9w 

w9r 

Tzz 

n9r  rJ 

* Tz0 

= 0. 
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To  calculate  numerically  the  stress  components,  the  analytical 
expressions  must  be  first  obtained  from  Equations  (23),  (36),  and 
(46).  The  calculations  involve  taking  derivatives  of  Bessel's  func- 
tions of  different  orders.  Thus,  the  expressions  for  the  stresses 
are  Bessel's  functions  of  various  combinations.  They  will  not  be 
given  here  but  will  be  entered  directly  into  the  Fortran  program 
written  for  the  determination  of  these  stresses. 


COMPUTATIONAL  RESULTS  AND  CONCLUSIONS 


Three  separate  Fortran  programs  are  written.  One  computes  the 
radial  displacement  u from  which  the  normal  stresses  Trr,  Tqq,  and 
tzz  are  computed,  the  second  program  computes  the  circumferential 
displacement  v from  which  the  shearing  stress  TrQ  is  computed,  and 
the  third  program  computes  the  axial  displacement  w and  the  shearing 
stress  Trz.  A sample  calculation  is  performed  for  a band  a = 2.034 
inch  and  b = 2.112  inch,  with  material  constants  E = 106  psi,  u = 

.26,  p = .000129.  This  sample  computation  is  based  on  an  angular 
velocity  history  with  u>q  = 60950  radians  per  second  and  tg  = 0.01 
second.  A second  calculation  is  made  for  a rotating  band  of  the 
same  inner  radius  but  larger  outer  radius  to  assess  the  effect  of 
the  increasing  thickness  of  the  band.  Results  were  plotted  to  show 
the  variations  of  the  stresses  with  the  radius  and  also  the  varia- 
tion of  stresses  with  time.  Figure  2 shows  the  three  normal  stresses 
plotted  vs  the  radial  distance  between  the  inner  and  the  outer  radii. 
It  is  evident  from  the  plot  that  the  most  severe  normal  stress  com- 
ponent is  by  far  the  radial  component  arr.  All  these  stresses  vary 
from  tension  at  the  inner  surface  of  the  rotating  band  to  compres- 
sion at  the  outer  surface  crossing  the  zero  value  at  the  midpoint  of 
the  thickness.  At  the  outer  surface  the  sample  calculation  gives  a 
compressive  stress  slightly  over  16,000  psi,  based  on  a material  with 
E = 106  psi  and  a geometry  corresponding  to  105  mm  gun.  A second 
sample  calculation  is  based  on  the  same  inner  radius  but  an  increased 
outside  radius  b = 2.2374  inch.  In  this  case  the  increase  of  radial 
pressure  is  very  considerable.  At  the  outer  surface  the  compressive 
stress  is  increased  to  about  90,000  psi. 

Figure  3 is  a plot  to  give  some  idea  about  the  increase  of  the 
radial  stresses  at  the  inner  and  the  outer  surfaces  of  the  rotating 
band  during  the  time  period  in  which  the  angular  velocity  of  the 
projectile  is  increased  from  zero  to  its  maximum  value.  It  can  be 
seen  that  the  increase  of  the  stress  is  less  rapid  than  the  linear 
variation  assumed  for  the  angular  velocity  vs  time  function.  This 
sample  variation  is  quite  typical  for  all  other  stress  components. 


•««««  „ d 

^:S:  ;i!£^  &> 

Vlz->  stress  ?“  ® °f  the  boundary  * fhtly  lower  than  ,1  h®  n°rmal 
components  co^itions  . £ fh°nditions  used t"?31  thoy 

--S'issiSH?  Swg 

To  refin,  th  ln"'r  a"d  °“‘*r 

hnn«  j_  ° tnOSP  . 


“cu  at  the  inner  anH  CCause 

To  refine  th  d °Uter 

boundary  condit  hCSe  results  it  io 

carried  out  1°^  be  assumed  in  recommended  that  non  ho 
because  thp  • much  less  Pff  ^ cases.  Thes^  „ n homogeneous 

Elated  in  the:gceralUeS  -ce-r^than  is  in^eS^J^0"*  ~n 

the  current  analysis*  the  computation have  ^ ana1^ 

tion/f'^^ion,  the  ..  ^ ^ 


-i-“oauon  have  b - 

In  conclusion  th 
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APPENDIX  A 


To  evaluate  the  normalization  constant  N defined  in  Equation  (15) 


N 


b ■ 

r*2(In  f)dr 

a 


(Al) 


it  is  recalled  that  ip  is  defined  by  Equation  (11);  thus  satisfies 
Equation  (6),  i.e.. 


dil  , I 4*  * 

dr2  r 37 


9 


which  can  be  recast  in  the  form 


r + Ce2'“T^  = °- 


(A2) 


Equation  (A2)  can  be  further  changed  into  the  form  (Ref  2) 


Therefore, 


and 


d_(rdi,2  , s2r2  dgil  . igl  , „ 


e2  fb  r2  , ((r^)V)b  , 0. 

J a 
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,IUBD 


Integrating  by  parts,  we  have 
b 


2e' 


,2  r diK  2 , 2 2 


r*  dr  = [r  (^r)  - (6  r‘-l)*]  . (A3) 


For  g- values  which  are  solutions  of  the  eigenvalue  equation  ip  (a)  = 
ijj(b)  = 0,  Equation  (A3)  is  reduced  to 


? 

rip  dr 


rW(g)2  - a2(g)2  ] 
s"  »u 


where 


Since 


26 


d^(6nr) 


1I[b2^’2(6nb)-a2^,2(6na)]  , (A4) 


'P  (Bnr) 


— ai^-  - 

Jl(6nr)  ' SnJ0(Bnr)‘7Jl(6nr)>  etc" 

6n  I W> 'KY0(6nr) > T?f J1 C8„r) (Snr) J . (A5) 


APPENDIX  B 


From  the  formulas  available  in  Reference  (3)  it  is  found  that 
b 


f «2>K6nC)d5  - ^-[52t2(Sn5)] 


b 

I 

a 


(Bl) 


where  <i'2(BnC)  = J2(Bn£)-K  Y2(BnC)  and  K is  the  same  constant  as 
appeared  in  Equation  (11). 


APPENDIX  C 


Eigenvalues 

The  eigenvalues  for  the  radial  displacement  problem  are  com- 
puted from  Equation  (10)  with  the  Fortran  program  Root  and  are  given 
in  the  following  table. 


k 

n = 1 

2 

3 

4 

5 

1.03834 

81.92742 

163.85156 

245.77344 

327.67969 

490.00000 

1.10000 

31.42676 

62.83725 

94.25140 

125.66644 

157.08179 

1.25000 

12.59004 

25.14465 

37.70706 

50.27145 

62.83662 

1.66667 

4.78508 

9.44837 

14.15300 

18.86148 

25.57148 

The  eigenvalue  equation  for  the  circumferential  displacement 
problem  as  given  by  Equation  (3)  can  be  reduced  to  the  following: 


I 


The  eigenvalue  equation  for  the  axial  displacement  problem 
after  some  algebraic  simplification  was  done  on  Equation  (45)  reduces 
to 

(kq)YQ  (q) -JQ  (q)Y1  (kq)  = 0 (C2) 


It  turns  out  that  both  Equation  (Cl)  and  (C2)  give  the  same  set 
of  eigenvalues  which  will  be  listed  below. 

k n = 1 2 3 4 5 

1.03834  40.5  123.0  205.5  287.0  369.5 


A 
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NORMAL  STRESSES 


n VI- Nr.  It:'  'MiTAI  15  > ,T(  20)  . R A ( AO  ) * SUMA  ( 151 1 SUKBl  151  ,SUKC( 
IT  All'll  10,  AH  I , T ALIT  I 10 , AO  1 , T AUZ  ( 10,  AO) 


0CC2 

HI:  All  I 5 , 2>! ) 1 1)  1;  T A 1 1 1 , 1*2,6) 

CC03 

20  rCRMATI 5T1C.51 

CCOa 

A = 2 • 0 2 A 

('005 

Rl‘0=  . OuO  1 29 

CCi'6 

SUVA  1 l 1 *0 

CC07 

SLMITI  11*0 

CCC8 

-SLVCI  11=0 

CC09 

CC1=?017190.0 

CC1C 

CC2  = A 30000 .0 

CC11 

CC3-1687300.0 

CC1J 

-RPS  = 60955-.XL 

CC  13 

TZ= . 0 1 

CC1A 

_C  = 2 . * I RPS/TZ  )-*.*2 

CC15 

CV=  1 . 038  3A 

CC  16 

.SCK=CM*CK 

CC  l 7 

T ( 1 1 =0 

ccls 

-R  A-C1-).=A_ 

cci9 

CC  50  1=2,11 

CC2C 

-T  1 I 1 =T  ( 1-L)  + »001_ 

CC21 

CC  50  J = 2 • 11 

CC22 

RAI J)=RA( J-ll+.OJ 

CC23 

CC  30  K = 2 , 6 

CC2A 

■ARl=aEIA.UC 

CC25 

6R2=CM*ARl 

CC26 

C=.0001 

CC  2 7 

CALL  HESJ(AR1,0,0J01,0,IERI 

CC20 

CALL  HESJ(AR2,0,BJ02,0,.IER) 

CC29 

CALL  BESJIAR1 . l.UJll.O, IER) 

CC30 

CALI — UESJ  LAR  1 ,X_,  BJZUO-.-l  ER  ) 

CC3 1 

CALL  BESJIAR2,2,BJ22,0,IER) 

0C32 

CALL— BESY.IARl.O, UVOl.O,. IER) 

CC33 

CALL  BE$Y(AU2»0,BY02»0»1ER) 

CC3A 

CALL  BESYI AR1, 1 , B Y 1 l , 0 , IER) 

CC35 

CALL  BESYIAR1,2,BY21,0,I6R) 

XC36. 

CAU — BEC-Y  I AR2.  2,  BY22.C.  I ER1 

CC37 

CT=BJ11/BY11 

CC38 

-X2=BJ21-CT*8Y2i 

CC39 

X2K=BJ22-CT*0Y22 

CCAC 

X1P1=BJ01-CT*BY01 

CCA1 

X1P1S=X1P1**2 

CC  A2. 

X IP  1 K = B J 0 B Y02 

CCA  3 

XIP1KS=X1PIK**2 

CCAA 

CU=SCM*X2K-X2 

CCA5 

AR1 A*AR 1/A 

CCA6 

ALN=(CG1/RHC)**.5*AR1A 

CCA7 

CC=ALN*»A*I SCM*X1PIKS-X1P1S) 

CCA8 

CN=C*0U/(Q0*AR1AJ. 

CCA9 

AA«ALN*TI  1 1 

CC5C 

BE* AA**2 • 

CC5  1 

CC-COSIAA) 

CC52 

T F «BB-2* 1 1 . *DD ) 

OC  53 

AR  3 = AR 1 *R A 1 J ) / A 

CC5A. 

-CALi—BE-SJ  I AK3..  0 ,.B  JOJ  *0  ._L£RJ 

CC55 

CALL  OESJI AR3, 1 , BJ13.0, IER  ) 

27 


) 


C C.  o 

CC*>  7 

C 0 ? ii 

CC  V7 

CC60 

CC61 

GC62 

CC63 

CC64 

CC  65 

CCfct 

CC67 

CC6S 

CC69 

CC70 

CC71 

CC?2 

CC73 

CC74- 

CC75 


call  bfcSYl  AR  J,0,h‘Y03,D,  IER  ) 

CALL  OESYt AR3, l,0Y13,O, IER) 
UPJ(I1J13-CT*0Y13I/RA(J) 

CUR=(BJ03-CT*»Y03l*ARlA 
FA=C01*DUK+C02*UR 
Ft2  = C01*UR*-CG2*DUR 
-FC*CG2*-(  OUR+UR-) 

CNTF=CN*TF 

■SLIKA  ( K ) *SUMA  ( K— -1-1  ♦CNTF*F  A 
SLPfll K )-SUMB(K-l )*CNTF*FB 
SCKC(K)  = SUf'C<K-l-)*CNTF*FC 
30  CCNTINUE 

-T-A  1>R  H-t- J4-»-$UMA-I  6-1 
T MJT ( 1 , J)=SUMB<6) 

TAUZII  ,J)-*SUMC<61 

KR I TE ( 6 , AO ) I * J , T AUR ( I * J ) * TAUT ( I , J ) , TAUZ ( I , J ) 

<•0  FCRMAT-I  lHtlXi-l.Z+-l3t5X+£lS*5r*XfEl5r5,SXtEl5.S)- 
50  CCNTINUE 
GALl-EX-d 
ENO 


I 


CIRCUMFERENTIAL  STRESS 


'ClU 
CCC2 
CCC3 
CC  O'. 
CCC:» 
CC06 

reel 

CC'-c 
CCC9 
CC10 
CCI1 
CC12 
CC1  3- 
CCH 
CC  1 5 
CCle 
CC  1 7 
CClfl 
CC  19 
CC2C 
CC2-1 
0C22 
CC  23 
CC2A 
CC25- 
0026 
CC27 
CC  23 
CC29 
CC3C 
CC31- 
CC32 
CC33 
CC3A 
CC36- 
CC36 
CC-37- 
OC  33 
CC  39 
CCAO 
CC41 
CC42 
CCA3- 
CCA4 
CCA5 
CCA6 
CC47- 
CCA8 
■ C C A9- 
CC3G 
CC51 
CC  52 
CC53 
0C54 
GC5S- 
CC  56 


CIMENSICN  PcTA(  15)  ,T120) , R A < 554 , SUMA < 154 , TAUZT ( 10 »60 > 
K E AO  ( ’> , 20 ) ( 0 E T A ( I ) • I * 2 • 0 ) 

20  rCRMATI7E6.nl 
A * 2 • 0 2 A 
m-Gs. 000129 
CC3»793650 
RPS-60955,0 
TZ* . 0 1 
C*2.*RPS/TZ 
CM* 1 .03836 
SCM*CM.*CM 
T( 11-0 
Wt-l-W 
sumai  i j =o 
DC  M 1*2,11 
T ( I >=T( l-l J+.OOl 
CC  50  J = 2,  11 
RA{J)=RA(J-l>f.01 
-OC— 30  -K*2,8 
AR1=BETA(K) 

-AR2=CM#AR1 

C=.000l 

•CALL— BES  J (AR1*0»8J01»D, IER) 

CALL  BESJ(ARljl,BJll,0,IER> 

-CALL— BE  SJ  (-AR2*  l t-BJ  1 2-rO-r-I  ER  1 
CALL  BESJI AR1*2.BJ21,D, IER) 

CALL  0ESJ(AR2*2»0J22,0,IER) 

CALL  BESY(AR1,0,8Y01,0,IER) 

CALL  RCSYURl,  l,BYll,D,IER) 

CALL  BESYI AR2, l, BY12.D, IER  ) 

CA  L-  L-— 0 E S Y(  A R l_,-2-rC-Y-2  1,0  IER4- 
CALL  BESY(AR2,2,BY22,0,IER) 

X2  = B.)21*BYU-BJU*BY2l 

X2K  = BJ2  2*BYll-tBJll*BY?7 

CU=SCM*X2K-X2 

#1P1=BJ01*BY11-BJ11*BY01 

•X  1 P 1-S*X-1  Pl**2-I 

Xlf>lK  = BJl2*8Yll-BJl  1*BY12 

X1PIKS*X1P1K**2 

AR1S=AR1*AR1 

CC=(SCM,-1/AR1S)*X1P1KS-X1P1< 

AR I A*  AR  1 / A 

ALN*7BA  37*AR  1 A 

CN=C*CU/(«C*ALN*AR1S) 

TF»1-CCS(ALN*T(  I ) ) 

AR3*AR1*RA< J)/A 

CALL  BESJ< AR3,0,BJ03,0,IER) 

CALL  HESJIAR3, 1 , B J 1 3 , 0 , IER) 

CALL  - UCSY  WR3.0,BY03fC«-ICR-) 

CALL  BESYI AR3, 1,BY13,D, IER) 

VR* ( B J 1 3~CT*BY13 ) /RA ( J ) 

0VR1*AR1*(BJ03-CT*BY03) 

CVR2  = - ( B J 1 3-CT*BY1.3  )-ARA(  J4 
0VR*0VRl+0VR2 
•RA-CC3RI OVU-VR 
CNTF»CN*TF 


29 


t 


CC'j7 

A C K ) af.Uf'A  (k-n+CNTF*FA 

ocr>e 

30 

CONTINUE 

CC5S 

TAUZTf  I » J ) *=  SUHA  ( 6 ) 

CCoC 

KKITE(6«A0)ItJ»TAUZT(IfJ) 

CC6  1 

AO 

FORMAT  ( IX , 1 3 ,!l  3 , 5X , 6 15 .53 

CC62 

50 

CONTINUE 

CC63 

^.ALL-eX-l-T 

CC64 

END 

AXIAL  STRESS 


CO  Cl 
CO  02 
CC03 

rco'i 

CCC5 

ccct 

-CG07 

ccoo 

cco<? 

CC10 

cell 

CC  1 < 

CC  1 3 

CC1A 

OC 1 5 

CO  16 

CC  17 

CC  1 8 

CC1.9F 

CC  20 

CC21 

CC22 

CC23 

CG2A 

0C-25 

CC26 

CC27- 

CC28 

CC29- 

CC30 

CC31- 

CC32 

CC33- 

CC3A 

CC35 

0036 

00  3 7- 

CC38 

CC39- 

CCAO 

OCA  1- 

CCA2 

OCA  3 

CCA  A 

CCA5 

CCA6 

CCA7 

CCA  8 

CCA9 

CC50 

CC51 

CC52 

CC53 

6C5A 

-OC-5-5 

C096 


«;  I VI- ,MS  I r-N  Hi  TA( 19) ,T120),RA(**),SUMA<19>tTAUM  10,  AO) 
KLA0I5,2U)  (BETA!  11,1-2,8) 

20  FORMAT  I 7F6.1  ) 

A - 2 • 02  A 
«FC=. 000129 
003=793690 
RPS=60955-«0 
TZ-.01 

0 = 2.  *RPS/. 13165*72 
CK  = 1 • 0303A 
SCM=CM*CM 
T ( 1 ) =0 
-R  AI-l  )-=A 
SUM A { 1 1=0 
CC  50  1=2,11 
T ( I ) = T( 1-1  (♦.001 
-CC  50  J = 2 « 1 1 
RAI  J)=RA(  J-D  + .Ol 
-CO-  30-  K=2,8 
ARl-OETA(K) 

AR  2 = CM* AR 1 
C=.0001 

CALL-OESJI AR l, 0,8 JO 1,0, IER) 

CALL  BESJ ( AR2,0,BJ02,D, IER) 

-O-ALL — llE  S-J  ( AR  1-r-l-rtW  l-l-r&.-I-CR-) 

CALL  BESY(ARl,0,BY01,0,IGR) 

-CALL  BESYt AR2,0,8Y02,D, IER) 

CALL  BESY(AR1,1,BYU,0,IER) 

CT=BJ01/BY01 
X2=BJ1 l-CT -0YO1 
-X2K=8  JO  1-CTABY02- 
CU=SCM*X2K-X2 
-X1P1-BJ1 1-CT*BY01 
X1P1S=X1P1**2 
X1P1K=DJ02-CT*OY02 
X1P1KS=X1P1K**2 
-Aft  1-S-—A  ft 1 —A  ft  1- 
CC=SCM*X1P1KS-X1P1S 
-AR  1 A=  Art  1 / A 
ALN  = 78A37*AR1  A 
CN=C*0U/(0D*ALN*AR1S) 

TF-1-C0S1 ALN-Tl  I ) ) 

-AR3=AR-1*RA(J  )/A 

CALL  BESJ! AR3,0,BJ03,0, IER) 

CALL  BESJ  (AR3,  1,BJ  13,0, -IER  ) 

CALL  BESYl AR3,0,BY03,0, IER) 

CALL  BESYl  AR3,  1, BY  13,0,- I-ER4 
CKR1=AR1A*(8J13-CT*BY13) 

-OfcR2«~<8J03-  C-TRB-Y43-)  V-R  A (-J-) 

ChR«-DWRl+0Wh2 

FA»C03*0WR 

CKTF«CN*TF 

SUMA ( K ) -SUMA  ( X.-*l  > *CNI£*FA 
30  CONTINUE 

T-AUZ-l-I-r-J)— -SUMA-t-ft-1 

WR  ITE  (6,  AO)  I ) J.TAtJZI  I . J) 


31 


r 1 


ccr>7 

AO 

FORMAT  ( lXf-I3t  l3»  5X,E15.5  J 

CC56 

50 

CONTINUE 

0C5S 

CALL-EXIT 

CCM'. 

ENO 

\ 


I 


CCi'  ’ 

;cr 1 

0CC4 
CCCs 
CCCd 
C C C 7 
crco 
CCCS 
CCIC 
CCll 
CC 12 
CC 1 3 

cca<t. 

CC  1 !> 

CC  16 

CC  1 7 

CC1Q 

CCI  9, 

CC20- 

CC21 

CC22. 

CC2  3 

CC24 

CC25 

CC26. 

CC27 

CC26. 

CC2S 

CC3C 

CC  3 1 

CC32. 

CC33 

C C 34. 

CC35 

CC36 

CC37 

cc3  a_ 

CC39 

0C4C. 

GC4  1 

CC42 

CC43 

CC44. 

CC4  ■> 

CC46 

CC4' 

CC-.fc. 

0C49 

CC50- 

CC  6 1 

CC52 

CC  Vi 

CC54 

CC55 

CCS  6. 


ROOT 

IC.N  X(  1000),  OX  C 100>»$UCOO).DF(lCOO)tCNMr.nfti 

X ( l ) =0. 

X(2)  = 4 1 . 

CF ( 1 ) = 1 .0 

rx (2)=.s 

AR1=X<2) 

AR2=1 .03834*AR1 
C*.0001 

CALL  ttESJIARl,  l.UJll.D,  IER) 

CALL  BESJI AR2,0,BJ02,0, ItRI- 
CALL  HESJI AR2,  1.BJ12.D. I ER ) 

CALL  BES  Y( AR 1 , l.BYll.O. ICR ) 

CALL  HESY(AR2.0.BY02.0.IER) 

CA  L L — B £S  Y.U  R i fC+4-EJU 

C = UJll/OYU 
FA«BJ02-C*BY02 
FC»BJ12-Q*BY12 
CFI2)»FA*AR2-FB- 
OC  15  J*1.1000 
•CMJ-Ufl- 
15  CONTINUE 
S (3 ) =0 

CC  10  1 = 3. 1000 

IFICNII).  EC.„1  .L-GO-TC  35 

IF(A8S(DF( 1-2)  ) . LT.  .00011  GO  TO  35 

XC.= X L L--CL=JCLLr_2J 

I F ( S ( I > . EC.  1.)  CO  TO  45 
IF(AOS(OF(.L-1)  ).-LT._.0001>  G0_ TO  35 
oxi  n=cx(  i-n 

GC  TO  65 

45  OX  ( I ) =0X  ( I-D*. 5 

A 4 XFXS  IJJ EIU_1— ) .GO  .TD..S5 

I F { XD . GT.  0. ) GO  TO  35 
XU.)  =X(-L=OX=DXLU 
GC  TO  20 

_35_XU  )-X(  I.-1J.*DX(IJ 
GC  TO  20 

-55..  UU.X  0 . GT  ...  0 . J_GQ— T 0—75 

XI  I ) = X< 1-11*0X1  I ) 

GC-TQ.  20 

75  X(I)-X(I-1 >-0X ( I ) 

GC  TO  20 

'20  AR1?X<  I ) 

AR2*i.Q3B34*ARl_ 

CALL  UESJIAR2.0.BJ02.0.IER) 

CALL-OESJI AR 1,1, BJ 11,0, IER) 

CALL  OESJI AR2, 1.BJ12.0.  IER) 

.CALL. BESY(AR2,0,BY02, D, IER) 

CALL  UESYIARlil.BYlltO.IER) 

CALI BES Y ( AR2 , 1 , B Y 1 2 , Q . LEJL1 

FA=BJ02*BY11-BJ11*BY02 
FB  = BJ12*BYU-BJ11*BY12 
OF  I I )»FA*AR2-FB 

.1 F ( ABS  ( OF  ( I )).  LT...000L1.  GO -TO- 100 
OC-OFI I)*DF( 1-1 ) ■ 

IFIOO.  GT. -0.1. GO. TO. 90. 
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